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DENSITY OF EIGENVALUES OF RANDOM NORMAL 

MATRICES 

PETER ELBAU AND GIOVANNI FELDER 



Abstract. The relation between random normal matrices and conformal map- 
pings discovered by Wiegmann and Zabrodin is made rigorous by restricting 
normal matrices to have spectrum in a bounded set. It is shown that for a 
suitable class of potentials the asymptotic density of eigenvalues is uniform 
with support in the interior domain of a simple smooth curve. 



1. Introduction 

In recent work initiated by P. Wiegmann and A. Zabrodin [1, 2, 3, 4], a connection 
between the normal matrix model and conformal mappings was discovered. In this 
model one considers random normal N x N complex matrices with probability 
measure 



CN 

>. 

TJ- . (1) P N (M)dM = Z„ 1 exp{-(N/to)tr(M*M-p(M)-p(M)*)}dM, 

o ■ 

\Q ■ where dill is a natural measure on the variety of normal matrices and Zjy is the 

*vO | normalization factor. The result is that, as N — > oo, the density of eigenvalues is 

1/nto times the characteristic function of a bounded domain in the complex plane. 
This domain is characterized by the fact that its exterior harmonic moments are 
the coefficients tj of the polynomial p appearing in the measure. Moreover, the 
Riemann mapping of the exterior of the unit disk onto the exterior of the domain 
obeys, as a function of the tj , the equations of the integrable dispersionless Toda 
hierarchy. 

These fascinating results remain however at the level of formal manipulations 
of undefined objects, as the integrals diverge, except in the simplest case of a 
S^ ■ polynomial p of degree 2, where the domain is bounded by an ellipse. 

?H ' The purpose of this note is to give a setting in which the above statements make 

mathematical sense and to give a proof of these statements. 

The problem of divergence of the integral over normal matrices is solved in a 
naive way, by restricting the integral to normal matrices whose eigenvalues lie in 
a compact domain D of the complex plane. Then for small to the results can be 
formulated in terms of polynomial curves, i.e., curves in the complex plane admitting 
a parametrization of the form w h- > h(w) = rw + $^--_ a j w ~ J i \ w \ = 1- For simple 
polynomial curves the problem of determining the curve out of its exterior harmonic 
moments tj = -^ § zz~i dz and the area irto of the interior domain (the interior 
domain is the bounded connected component of the complement of the curve) has 
a unique solution for small to, as we show in Section 5. Our main result is then: 

Theorem 1.1. Let D C C be the closure of a bounded open set containing the 
origin. Let p(z) = tiz 2 + ■ ■ ■ + t n +\z n+l be a polynomial such that \z\ 2 —p(z) —p(z) 

l 
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has a non- degenerate absolute minimum in D at z = 0. Then there exists a 5 > 
so that for all < to < 5, 

(i) there exists a unique simple polynomial curve 7 with exterior harmonic 

moments t\ = 0, £2, • • • , t n +i, 0,0,... and area of interior domain irto; 
(ii) the expectation value of the density of eigenvalues of random normal ma- 
trices with spectrum in D and distribution (1) converges as N — > 00 to a 
uniform distribution with support in the interior domain of 7. 

The condition on p implies the Hessian condition \t2\ < 2 an d it is fulfilled if 
the Hessian condition holds and tz,...,t n are sufficiently small. It then follows 
from results of [5, 1] that the curve 7 as a function of the tj in this range provides 
a solution of the intcgrable dispcrsionless 2D Toda hierarchy obeying the string 
equation. 

The paper is organized as follows: the basic definitions of the random normal 
matrix model are recalled in Section 2. We then introduce the "equilibrium mea- 
sure" as the unique solution of a variational problem in Section 3 and show in 
Section 4 that the density of eigenvalues converges to it. These are either known 
results or adaptations of results known for hermitian matrices to our case. In Sec- 
tion 5 we introduce the notion of polynomial curve and prove Theorem 5.3, which 
is a stronger form of part (i) of Theorem 1.1. In Section 6 we prove part (ii) of 
Theorem 1.1, sec Theorem 6.1, and discuss our results in Section 7. 

2. Eigenvalues of random normal matrices 
We consider the probability measure 

P N {M) dM = -L c -NtrV(M) dM) Z N = [ P N {M) dM, 

defined by a potential V, on the set 

Af N (D) = {M e Matc(iV) I [M,M*] = 0, a(M) c D} 

of normal N x N complex matrices with spectrum in some compact domain DcC. 
The measure dM is the Riemannian volume form on (the smooth part of) AFn(D) 
with respect to the metric induced from the standard metric on the vector space 
C N of all N x N matrices. In a parametrization by eigenvalues and unitary matrices 
it is given by [6] 

N 

dM= dU Y[ \z i -z j \ 2 l[d 2 z i , 

l<i<j<N i=\ 

where M — C/diag ((z^fL-^U* and dU denotes the normalized U(iV) invariant 
measure on XJ(N)/XJ(1) N . 

This leads to the probability measure 

N 1 N 



i=l l<i<j<N i=\ 



N 

2. 



■J D 1 ^ : ^ • ^- jv t a 1 



Zi 
l<i<j<N i=l 

on the space of eigenvalues z, e D, 1 < i < N. 
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3. The equilibrium measure 

We are interested in the behavior of the function Pn{z) as TV — * oo. Because 
the probability that two eigenvalues are equal is always zero, we may consider P/v 
as function on the set Dq = {z E D n | Zi ^ Zj\/i ^ j}. 

Introducing the probability measure 

1 N 



on D (xa shall denote the characteristic function of the set A), we write for z G D^ 
1 



N 



Pn(z) = —cx P \-N'U V(()d5 z (() + JJ \og\C-^\- L d6 z (Od8 z (0 

Letting N — > oo, only the infimum of the coefficient of —N 2 , 
Io := ^^ J ^)' 

1(H) := f V(z) dft(z) + ff log \z - CI" 1 d M (C) dfx(z), 
J JJz^c 

will be relevant. Here A4(D) denotes the set of all Borel probability measures on D 
without point masses (a measure with point masses could only arise from measures 
S z with some z% — Zj, but then P/v(z) = 0). Therefore, we safely can neglect the 
restriction on the double integral. 

The precise sense in which the infimum of I controls the large N behavior of 
Pjv will be discussed in the next section. Here we consider this reasoning only as a 
motivation for introducing the variational problem. 

Because of 7(&ftA) < oo (A the Lebesgue measure on C and \D\ = \{D) the area 
of D), Io is finite. 

Definition. An equilibrium measure for V on D C C is a Borel probability measure 
\i on D without point masses so that I(fi) = Jo- 
Theorem 3.1. Every continuous function V on a compact subset D has a unique 
equilibrium measure. 

In the remaining of this section we prove this theorem, which is a known fact from 
potential theory, see e.g. [7], and give necessary and sufficient conditions for jj, to 
be an equilibrium measure. The constructions are adapted from the corresponding 
results for Hermitian matrices, see [9]. 

3.1. Existence. To show the infimum is achieved, we choose a sequence (fJ, n )^Li 
in M(D) with I(p n ) — ► Jo- 
Lemma 3.2. The space of all Borel probability measures on D is sequentially com- 
pact. 

Proof: By the theorem of Riesz-Markov each Borel measure /jonD corresponds to 
exactly one positive, linear functional </> M € C(D)* and by the theorem of Alaoglu, 
the closed unit-sphere in C(D)* is weak-*-compact. Therefore, for each sequence 
(Mn)5£Li C A4(D), the sequence (<^„)5£Li contains a weak-*-convergent subse- 
quence (0 MnW )!£=!, i.e. 

3 </> G C(D)* : <^„ (fc) (/) - 4>(f) V/ G C{D). 
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Now we find a measure /jonD with (f> — (p^. This measure fulfills 
fd» n{k) ^ffdn (fc-oo) V/eC(£>), 

and hence is again a Borel probability measure. □ 

Because of this lemma there exists a convergent subsequence (n n tk)) of (fJ, n ) and 
a Borel probability measure fj, with /i rl (fe) — > M- 

To prove that /(/i) = Iq; we estimate with an arbitrary real constant L 



lim I{/J, n (k)) = Ihn / V(z) d/w fe )(;s) + lim \og\z - (\ d^„ (fc) (C) d/^ n ( fe ) {z) 

k— »OC k—>QO J k— »CO J J 

> / F(z) d^(z) + lim // minjlog |z - C| _1 , L} d / u„( fe )(C) d/x n(fc )(^). 

Approximating uniformly the second integrand according to the theorem of Stonc- 
Weierstrafi up to e > with a polynomial in z, z, Q and C, and using Fubini's 
theorem, we get a lower bound for the limit: 

lim I{n n{k) ) > I V(z) dft(z) + fj minjlog \z - (\~ l ,L} dju(C) M?) - 2e. 

Letting first e — > and then L — > oo shows that ^ has no point masses (otherwise 
the right hand side would diverge) and Iq = I(fJ-). 

3.2. Uniqueness. Next, we want to show that there is exactly one measure /x € 
M(D) with I(p) = Iq. So suppose ji e M(D) also fulfills I (fl) = I . Then we 
consider the family 

fit = tfi + (1 - i)/U = /i + £(/} - /i), £ e [0, 1], 

in.M(Z>). 

Lemma 3.3. Let fi and p, be probability measures such that the function log |z— C| _1 
is integrable with respect to fi® /z and io /i(g>/2. Then log |z — C| _1 is o/so integrable 
with respect to fi ® /L 

Additionally, we have the inequality 

(2) || log |z - CI" * d(A - M)(C) d(/i - fi){z) > 0, 

wii/i equality if and only if fi = /i. 

Proof: We start with the distributional identity 

log \z\- l A<p(z) d 2 z = -2tt^(0) 

for any Schwarz function tp. Introducing the Fourier transform ip of ip, we find 
i-1a,./,\j2 / 1 ni„i2~n„\\ j2, 



log|z|- 1 A^(z)d^= - / — (|£;|^(A;))d^ 



|fc| 2 

^JJw ( c ^ {kz+ks) - f{k) ) d2fc A ^ (z) d2z ' 
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where / denotes a real, continuous function which is one in the vicinity of zero and 
becomes zero at infinity. So we have for all z G C \ {0} the equation 

log \z\~ x = ^Jj^I (e*<** +£2 > - /(*)) d 2 k + C(f) 

with some real constant C(f). 

So, with Tonclli's theorem, we see that 

2 



log Iz-Cr 1 d(£-A*)(C) d(ji-fi)(z) = ^ / ^|2 



8'< fc * +fc *>d(£ -/*)(«) 



A 2 k 



is non-negative and finite, which immediately implies the intcgrability of log \z — w\ 
with respect to fi ® /}. To achieve equality in (2), we need 

e3(**+ fe )d/K(*)= ( ^ {kz+kz) dfi{z) 



for all k £ C, which reads /i = jx. D 

So we can expand I(nt) and obtain 

I(jH) = I (p) + t J (V(z) + 2 J log \z - CI" 1 d^C)) d(A - M )(«) 
+ t 2 /'/log|^-Cr 1 d(M-M)(C)d(M-/»)W. 



Lemma 3.3 now states that the coefficient of t 2 is non-negative and so the function 
t i— ► I (fit) is convex on [0, 1]. In particular, 

I(l H )<tI(fi) + (l-t)I(n)=I , 

which implies I(p,t) = lo- This requires the last summand in (3) to vanish and so, 
again with Lemma 3.3, we see that fi = JX. 

3.3. A variational form. To determine if a given measure /j, is the equilibrium 
measure for the potential V on the domain D, we may check a variational principle. 

Proposition 3.4. The probability measure /i is the equilibrium measure for the 
potential V on the domain D if and only if the function 



(4) E{z)=V{z) + 2l\og\z-t\- 1 MQ 
fulfills the relation 

(5) I E{z)djx{z)> E(z) dii{z) =: E for all ft e M(D). 

Additionally, we have the property 

(6) E{z) = E /i-almost everywhere. 

Proof: Let us first assume /i is the equilibrium measure. Then, for an arbitrary 
measure Jx € A4(D), the condition 

^/(m*)L =0 >o, nt = tLi + (i-t)n, te[o,i], 

has to hold. This means 

(v(z) + 2 J log \z - Cr 1 d^C)) d(ji - n)(z) > 0, 
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which immediately implies equation (5). 

Assume on the other side jjl fulfills condition (5). Then we obtain with the 
equilibrium measure no 

/(Mo) = Hji) + J (V(z) + 2 J log \z - CI" 1 dM(C)) d (Po - H){z) 

+ jj log \ Z - cr 1 d( Mo - m)(o d (^o - n)(z) 
> Up), 

and so \i = n$. 

To show the additional statement, we consider for the probability measure (X the 

set 



B:={zeD\ E{z) < E }. 

71(B) 1 



If n(B) > 0, the variational principle (5) for the measure fi = ^, jj would yield 



E < j ' E{z)*^Mz)<E , 

and so n(B) has to vanish. Thus, we get condition (6). □ 

Instead of verifying condition (5), we prefer to use the following statement: 

Corollary 3.5. If for a measure fi G M.{D) the function E, defined by equation (4), 
fulfills, for some real constant Eq, E(z) = Eq on the support of fi and E(z) > Eq 
everywhere, then jj, is the equilibrium measure. 

4. The eigenvalue density 
Definition. The fc-point correlation function R^ is given by 

TV! r N 

<\^)U) = 7W^yJ DN _ k M^)^) U d 9 *. 

So, the one-point correlation function is up to normalization nothing but the density 
of the eigenvalues. 

As indicated in the previous section, all the correlation functions can be cal- 
culated in the limit N — > oo out of the equilibrium measure /i, as in the case of 
hcrmitian matrix models, see [8, 9]. 

Theorem 4.1. For all <p € C(D k ) we have the equality 

k k 

(7) a^oo / -L^ z ^) r( n ((*)*=i) n <*** = / ^((«i)*=o n d ^- 

^°°JD k i=1 J i=1 

I.e. the measure jt^R n [(zi) k =1 ) Y\i=i ^ z i on D k converges weakly toYii—i d/j(zi). 

Proof: Substituting in the left-hand-side of equation (7) the definition of the cor- 
relation functions and turning our attention on the highest order in N, we obtain 
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(because Pjv is invariant under the symmetric group) 

iV TV 

i 1 ,...,i k =i jDk i=i 

Since for large values of TV the probability distribution localizes at values z <E L>^ 
with I(i5 z ) « lo, let us consider the sets 

^V„ = {* e £>^ I /(<5 Z ) < h + vh V> o. 

Lemma 4.2. 27ie probability P^(D N \ Aj^^) drops for N — > oo exponentially to 
zero. 

Proof: For an absolutely continuous equilibrium measure d/i(z) = V'C 2 ) d 2 z, we 
get with Jensen's theorem and J I{S Z ) f\ d[i(zi) = I + o(l) 

r N 

Z N > / e^ 27 ^)^™ i lo ^^) TT d M (z,) > c -^ 2 ^+o(^) ; 

and therefore, 

- N 

P N (D N \A Nr ,)< / e A f2 /o+o(iV 2 )-^(/ 0+ ,) TT d 2 Zt = o{c -N 2 v/2y 

Let the continuous function -jb ^2 0(( z ij )j=i) take its maximum on the compact 
set An. v at C, and set vjq,. q = 5^. Then, 

1 \ 1 N f k 

w R ») - w £ *«k)*=i) - / ^((^)ti) n 4"»m- 

ii.....ik — l ' i— 1 

Because of Lemma 3.2, we find a convergent subsequence i>N(n),-q ~ * v r) {n ~ * °°) 
with 

So (^ a^) - / ^w n dv ^ zi) - 

Lemma 4.3. We ftai/e ^ G A4(D) and, in the limit 77 — s- 7 I{v 71 ) — ► Iq. 
Proof: Using £ £ A N ( n y v , we obtain with the cut-off Lei: 
1 N{n) 1 _ 

V ; »=1 V ; l<i^j<JV(n) 



= / V(z)dl/j\r( n ),,,(z) 

+ / / minjlog \z - CI -1 , i} dv N(nhrl {() dv N(nhrl (z) - — — 



This case suffices for our needs, but the restriction is in fact not necessary. Indeed, we could 
perform an analogous argument for the measures dfi £ (z) = ip £ (z)d 2 z, ipe{z) = — % f B , > dfi. In 

the limit e — ► 0, where I(fi £ ) — > /oi this would yield the desired statement. 
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Sending first n and then L to infinity brings us to v 71 G A4(D) and therefore 
Jo < J(i/„) < /o + V- ' □ 

So, letting 77 — ► 0, a subsequence of v„ converges to the equilibrium measure /i, 
and thus, 

k 

jE^ U j^R^) < j 4>{{*i)U) II M*i)- 

* ' ' i— 1 

Arguing on the same way for the limes inferior concludes the proof. □ 

5. Polynomial curves 

Definition. A polynomial curve of degree n is a smooth simple closed curve in the 
complex plane with a parametrization h : S 1 C C — > C of the form 

(8) h(w) — rw + ao + a\w~ + • • • + a n w~ n , \w\ = 1, 

with r > and a n ^ 0. The standard (counterclockwise) orientation of the circle 
induces an orientation on the curve. We say that a polynomial curve is positively 
oriented if this orientation is counterclockwise, i.e., if the tangent vector to the 
curve makes one full turn in counterclockwise direction as we go around the unit 
circle. 

Proposition 5.1. Let 7 be a positively oriented polynomial curve with parametriza- 
tion h of the form (8). Then h, viewed as homolorphic map onC x , restricts to a 
biholomorphic map from the exterior of the unit disk onto the exterior 0/7. 

Proof: We have to show that h'(w) 7^ for all w in the complement of the unit disk. 
Let t denote the tangent vector map w 1— » t(w) = h'{w)iw = \{rw — ^jdjW^ ). 
Since 7 is a simple closed curve, the map w f— > t(w)/\t(w)\ is a map of degree 1 
from the unit circle to itself. Therefore we have 

1 = — i dargOH) 

27T J\w\=l 

1 r t'{w) 



r . dw 

2?ri / h=1 t{w) 

1 / ?M , 



27ri J\ W \ =R t(w) 

Here N > denotes the number of zeros of t(w), counted with multiplicity, in the 
complement of the unit disk and R is so large that it contains them all. The latter 
integral is 1 as can be seen by sending R to infinity. Thus ./V = 0, and h! has no 
zeros in the complement of the unit disk. □ 

A simple consequence of this proposition is that a polynomial curve is uniquely 
parametrized by a map of the form (8) with r > 0. Indeed, any other conformal 
mapping of the complement differs by an automorphism of the complement of the 
unit disk. But non-trivial automorphisms are given by fractional linear transfor- 
mations which do not preserve the conditions. 

From now on, we will only consider polynomial curves encircling the origin, i.e., 
such that the origin is contained in their interior domain. This can always be 
achieved by a translation, i.e., a shift of clq. 
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Definition. The harmonic moments (tj)j?-i of the exterior domain D- of a polyno- 
mial curve (or more generally of an analytic curve) encircling the origin are defined 

by 

£,• = / z~' J A 2 z = — — — <z> zz 



3 nj J D 2ttij j 1 

where only the right integral should be taken as a definition for j '• < 2. 

Proposition 5.2. Let 7 be a positively oriented polynomial curve of degree n en- 
circling the origin. 

(i) The exterior harmonic moments tj of 7 vanish for all j > n + 1 . 

(ii) There exist universal polynomials Pjj~ G 7i[r,aQ, . . . , dk-j], 1 < j < k, so 
that for j = 1, . . . , n + 1, 

n 

(9) jtj = a,j-ir~ j+1 + J^ akr~ k Pj,k(r, a , . . . , dk-j)- 

k=j 

Moreover, Pj.k is a homogeneous polynomial of degree k—j + l and it is also 
weighted homogeneous of degree k— j + 1 for the assignment deg(aj) = j+1, 
deg(r) = 0. 
(iii) The area of the domain enclosed by 7 is irto where 

n 

t = r 2 -J2jWj\ 2 - 

Proof: Since 7 encircles the origin, h(w) never vanishes for \w\ > 1. Hence the 
contour in the formula for tj may be computed by taking residues at infinity. For 

jtj = ■— (f> h{w~ 1 )h'{w)h(w)^ : ' dw 




The integrals in this sum vanish if k < j — 2. The formula for tj in terms of ak, r is 
obtained by expanding the geometric series and picking the coefficient of w _1 in the 
integrand. This proves (i) and the first part of (ii). The homogeneity property is 
clear. The weighted homogeneity follows by rescaling w in the integral. The same 
formula can be used to compute to, but the first term no -1 in /i(w _1 ), which does 
not contribute to the integral and was omitted for j > 1, must be added here. □ 

Examples. The terms in tj involving polynomials Pj t k with k < 3 are 

1 a 2 (2 air — a 2 ) + r~ 3 a 3 (3 a air — 3 a 2 r 2 — a 3 ) + • • 

3 r~ 3 a 3 (air — a 2 ) H 



Theorem 5.3. Let £2, • • • , tn+i be complex numbers so that \t2\ < 1/2. Then there 
exists an Aq = Ao(t2, ■ ■ ■ , £«+i) > so that for all A, t\ with < A < Aq and 
\ti\ 2 < A(l/2 — I £2 1), there exists a unique positively oriented polynomial curve 
of degree < n encircling the origin, with area A and exterior harmonic moments 
ti, ■ ■ ■ , £ n +i,0, 0, . . . 



h -- 


= a - r~ 


_1 aiao — r~ 2 


2£ 2 = 


= r~ a\ - 


- 2r~ 2 a 2 a 


3£ 3 = 


= r~ 2 a 2 


- 3r _3 a 3 a 
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Proof: The idea is to invert the map (r, ao, • • • , a n ) l— * (to, • ■ • , t n +i) for small r and 
ao- Set aj = r~ 3 aj, p = r 2 and consider instead the polynomial map 

F : (p, ao, • • • , a„) — » (to, • • • , tn+i), 

as a map from M x C" +1 to itself. The first claim is that this map has a smooth 
inverse in some neighborhood of any point (eRx C n+1 such that to = t\ = and 
|<2 1 7^ 1/2- By Prop. 5.2, this map is given by 

n 

to = p-^fJJWA 2 

n 

J*j = ctj-i +y^^akPj,k(r,ao,rai, . . . ,r k ^ 3 ak- ) 
k=j 

n 

= a.j-i + 22 &kPj,k{p, ao, ai, . . . , otk-j)- 
k=j 

From the contour integral representation of the harmonic moments we get the 
integral 

Pj,k(p,a ,. .. ,otk-j) 

= — I w k -3 ( 1- f^lai^w- 1 - 1 | f 1 + Va,/^-'- 1 I dw 

for any sufficiently large R. By computing the residue at infinity, we can calculate 
Pj t k and thus tj up to terms of at least second order in ao, p, 

to = p(i-K| 2 ) + --- , 

jtj = 6tj-i - j Ojao - (j + l)pa j+ iai -\ , j > 1. 

Hence F(0, 0, 2 £2, ■ ■ ■ , (n + 1) £™+i) = (0, 0, t2, . . . , t n +i) and the tangent map at 
this point sends (p, do, ■ • • , d n ) to (to, ... , t„) with 

to = (l-4|t 2 | 2 )p, 

jij = «j-i - it? + 1) tj+ido - 2j(j + l)(j + 2) t ]+2 t 2 p, j > 1. 

The tangent map is invertible if |t2| 7^ 1/2. By the inverse function theorem, F 
has a smooth inverse on some neighborhood of (0, 0, £2, ■ • ■ , t n +i). If (£2! < 1/2, F 
preserves the positivity of the first coordinate. 

In terms of the original variables, this means that given any t = (to, £1, £2, • • • , £n+i) 
with small to > 0, t\ and such that |£2| =/= 1/2, there is a curve w 1— ► h(w) with 
h(w) = rw + ao + ra\w~ l + • • • + r n a n w~ n and aj ~ (j + l)tj+i. It remains 
to show that if r > is small enough, h parametrizes a positively oriented simple 
closed curve containing the origin We first show that h is an immersion. Since 
h'(w) — r — ra\w~ 2 + 0(r 2 ) and lim r ^o a i = 2 £2, we see that as long as fal 7^ 1/2, 
h'(w) does not vanish on the unit circle. Similarly, we show that h : S 1 — > C is 
injective: we have 

\h(w) - h{w')\ 2 = r\w-w' + 2£ 2 («T 1 - w'- 1 )\ + 0{r 2 ) 
= r\w - w' + 2t 2 (w - w')\ + 0(r 2 ) 

But the expression in the absolute value can only vanish for w ^ w' if |2£~2| = 1 
which is excluded by the hypothesis. Moreover h{w) = rw + t\ + 2rt 2 w~ 1 + 0(r 2 ). 
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Therefore h parametrizes a perturbation of an ellipse centered at t\ . The condition 
on t\ is a sufficient condition for this ellipse to contain the origin. □ 

Example. Let k > 3 and let us consider curves with tj = for all j =/= k. Then 
dj = for all j 7^ k — 1, so that h(w) = rw^ 1 + auto' -1 . The relation between 
(t ,t k ) and (r,afe_i) is t = r 2 - (k - l)|a fe _i| 2 , kt k = Ofe_ir _fe+1 . This map is 
a diffcomorphism from the region < r < (fc — 1) |ofc_i|, which the condition for 
h(w) to be an embedding of the unit circle, onto the region 

< h < (k(k - l)|i fe |)- 2 /( fe - 2 )(fc - 2)/(fc - 1). 

As to approaches the upper bound for given t k , the curve develops cusp singularities. 

6. The equilibrium measure for a polynomial curve 
In this section we evaluate the equilibrium measure corresponding to potentials 

(10) V(z) = ±(\z\ 2 -2RcJ2tkZ k \. 

We anticipate the result: 

Theorem 6.1. For any set (tk)kLi C C with t\ = 0, \t2\ < \ and t k = for 

k > n + 1 and any compact domain D C C containing the origin as interior point 
and such that t$V , V given by equation (10), is positive on D\ {0}, there exists a 
5 > so that for all < to < S the equilibrium measure ji for V on D is given by 

1 

Trio' 

where D + denotes the interior domain of the polynomial curve 7 defined by the 
harmonic moments (tfc)£L > an< ^ ^ * s the Lebesgue measure on C. 
The rest of the section is dedicated to the proof of this theorem. 

6.1. The Schwarz reflection. We first need the notion of a reflection on an an- 
alytic curve (see e.g. [10]). 

Definition. The Schwarz function of an analytic curve 7 is defined as the analytic 
continuation (in a neighbourhood of the curve) of the function S(z) — z on 7. The 
Schwarz reflection p for the analytic curve in this domain is the anti-holomorphic 

map p{z) = S(z). 

Definition. Under the critical radius R of the polynomial curve defined by the 
parametrization h we understand the value 

R = max{|w| | h'(w) = 0, w e C}, 

which, by definition of the map h, is less than 1. 

Lemma 6.2. Let 7 be a polynomial curve parametrized by h and R its critical 
radius. Then the Schwarz function S and the Schwarz reflection p of the curve 7 
restricted to h(B 1 i R \ Br), where Br denotes the open disk with radius R around 
zero, are biholomorphic respectively anti-biholomorphic maps. They are given by 

(11) 5( ,) = *(_y and «z) = h(~±^ 

Therefore, p maps 7 identically on itself, h(Bi\Bn) to h(Bi/f>\Bi) and vice versa. 
Also, we have p 2 = id. 



M= 3-Xd+A, 



E(z) = V(z) + 4" / log 



z 

r 1 
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Proof: By definition of the Schwarz function, in a neighborhood of \w\ = 1, 

S(h(w)) = hiw- 1 ). 
Because for w <E <C\ Br the function h is biholomorphic, we may write 

^ = K/^))' ^ )= ^^G4), 
Taking the derivatives, we find that they do not vanish for R < |/i _1 (z)| < j^. D 

Lemma 6.3. In the interior domain D + of the polynomial curve defined by the 
parameters t k from (10) as its harmonic moments, the function 

is equal to zero and in the exterior domain D \ D + , its gradient reads 

(12) d s E(z) = Uz-p(z)). 

to 

Proof: To verify the first statement, we use Green's theorem and obtain 

\f iog|z-cr 1 d 2 c--ki 2 +Rc^ I ^io g |z-cr l c+^)dc. 

Integrating by parts of the second integrand yields immediately 

\ f log \z - CI" 1 d 2 C - -\z\ 2 - 2 Re ^- f log(C - z)(d(, 

and expanding the logarithm around z = leads us to E(z) = in D + . 

For the proof of the second part we write S = 5; + S c , where Si is analytic in 
D + and S c in the complement D \ D + . For the exterior function S e one finds with 
the Cauchy integral and Stokes formula: 

i fC-SiiC) A ,_ if i a2 , 



S c (z) = -—<b * >^' d( = - / -d% zeD\D + . 

27ri Jj C - z 7T ./£>+ z - C 

Because we know E to be constant on D + , 



1 / n+l 

d z E{z)^-[z~y j kt k z k - 1 -- — 

to\ fr[ *J D+ z- 



2, 



o = a ZJ e(z) = -!«->, *t fc «"- J - - / — — d^c I 

for all z G _D + . So on ini7, and by analytic continuation in the entire domain 
where S is holomorphic (which includes the exterior domain D \ -D+), we have 

Si(z) = E n ktlkt k z k - 1 . 

Therefore, for all z £ D \ D + , 

d- z E{z) = r( z - ^M - W) = T (* - P( z )) ■ 

t \ 'to □ 

Let us remark that this proof shows that the Schwarz function S has, at least 
around infinity, the form 

S(z) = J2 kt k z k - 1 + | + J2 v k z- k -\ 
fe=i fe=i 
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where the Vk = — J D z k d 2 z denote the harmonic moments of the interior domain 
D+. This fact was already used in [2] to establish a connection between the har- 
monic moments and the coefficients of the parametrization of 7. And as was shown 
in [2], we find, with the Theorems 4.1 and 6.1, that the Vk are nothing but the ex- 
pectation values of ^ tr(M fe ) with respect to the probability measure Pn(M) dM in 
the limit N — > 00 and hence, are completely determined by the harmonic moments 
of the exterior domain and to • 

6.2. The Gaussian case. In this case, where the polynomial curve is an ellipse, 
we are able to calculate the equilibrium measure on C explicitly. 

Proposition 6.4. The equilibrium measure \x for the potential 
V(z) = ±(\z\ 2 -t 2 z 2 ~i 2 z 2 ), |fe|<±, 



1 

ah 
where D + denotes the interior of the ellipse 



M = — Xd+\ 



,,,x Re(y^z) 2 Im(VM) 2 /1 + 2N, , / l-2|t 2 | JL 

Proof: As polynomial curve, the ellipse (13) has the parametrization 
(14) h(w) =r(w + 2t 2 w~ 1 ), r=- . 

We check that the given measure \x is the equilibrium measure by verifying the 
conditions of Corollary 3.5. To this end, let us introduce for \w\ > 1 the function 

1 d 2 C. 



£ (w) = E(h(w)), E(z) = V(z) + -^ r I log 

3+ 



*"*0 JD A 



z 

r 1 



Integrating equation (12) and its complex conjugate analog we get for £(w) the 
expression 

1 



£{w) = — I |/i(w)| 2 -|/i(l)| 2 -2Re / h{w- 1 )h'(w)dw 

Substituting in this expression relation (14) for h, we obtain 
t £(w) = r 2 (\w + 2t 2 w- l \ 2 - |1 + 2t 2 \ 2 ) 

- 2r 2 Re(t 2 w 2 + (1 - 4|£ 2 | 2 ) logw + t 2 vj- 2 - (t 2 + t 2 )) 
= r 2 (\w\ 2 - 1 - 4\t 2 \ 2 + 4|i 2 | 2 M- 2 + (1 - 4|t 2 | 2 ) log(| W |- 2 )) 

+ 2r 2 Rc(2t 2 ww~ l - t 2 (ww~ 1 )|u;p 2 - t 2 {ww~ 1 )\w\ 2 ) 
= r 2 (\w\ 2 - 1)(1 - 4|t 2 | 2 | W |- 2 ) + r 2 (l - 4|t 2 | 2 )log(| W |- 2 ) 

- 2r 2 (|w| 2 - 1)(1 - \w\- 2 )Rc(t 2 ww- 1 ) 

= r 2 (\w\ 2 - 1)(1 - 2|t 2 |)(l + 2|t 2 |H- 2 ) + r 2 (l - 4|t 2 | 2 ) logCM" 2 ) 
+ 2r 2 (M 2 - 1)(1 - |t«r 2 )(|*2| - Rc^ww- 1 )). 
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We are now ready to start estimating £{w) for \w\ > 1 and \t2\ < \- The last 
bracket we can estimate by 

\h\ - Refeww' 1 ) > \ti\ - \t 2 ww- 1 \ = 0. 
It therefore remains to show 

{\w\ 2 - 1)(1 + 2\t 2 \\w\- 2 ) + (1 + 2|t 2 | 2 ) logd^r 2 ) > 0, 
which follows immediately out of the following lemma. 
Lemma 6.5. For < a < 1 the function 

f(x) = (x- 1)(1 +ai _1 ) - (1 + a) log x 
is non-negative on the interval [l,oo). 
Proof: For the function / and its derivatives /' and /" we have 

/(1) = 0, /'(1) = and f"( x ) = ±(l + a-—)>0. 

So we showed that E(z) > for all z £ C \ D + . Because of Lemma 6.3, we 
also know that in the interior domain D + , E is zero and we therefore can apply 
Corollary 3.5 to see that /u is indeed the equilibrium measure for the potential V 

onC. □ 

6.3. The proof of theorem 6.1. As in the Gaussian case we are going to show 
that the measure /j, given in the theorem fulfills the conditions of Corollary 3.5 and 
is therefore the uniquely defined equilibrium measure. 

Because Theorem 6.1 is only valid for interior domains with small area, we are 
going to consider the asymptotical behaviour to — > 0, where the harmonic moments 
(tk)kLi are kept fixed. To catch the asymptotical behavior of the corresponding 
polynomial curve, let us parametrize it as in the proof of Theorem 5.3: 

n 

h(w) = rw + \ r J ajW~ J . 
j=o 

Then, for r — > 0, we have r 2 ~ to, ctj ~ (j + l)ij+i, j > 1, and, because we set 
t\ = 0, ao ~ r 2 . 

Lemma 6.6. The critical radius of h is asymptotically constant for r — ► 0: 

R= y/\aT\ +0(r). 

Proof: The roots of the function h'(w) = r — ra\w~ 2 — ... — nr n a n w^ n ^ 1 are in 
zeroth order at ±y/a± and [n — 1-times degenerated) at zero. D 

We consider now for z £ D respectively for w £ h~ l (D \ D + ) the functions 



(15) E(z) - V(z) + / log 

1 



-l 

j2 



cTC and 



/w 
hiw-^h' (w) dw 

We already showed in Lemma 6.3 that E(z) = in D + , so the first condition of 
the corollary is satisfied (this is essentially the way we have chosen our potential 
V). 
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Now, also with Lemma 6.3, we see that E > in the vicinity of the curve 7, 
strictly speaking in the domain h{Bi/n \ B\). Indeed, if we look at the connected 
components of the contour lines of the function E (which are smooth curves in the 
considered domain because there dgE 7^ 0), we see that the gradient vector d g E 
always points outwards, i.e. in the exterior domain of the contour line. Therefore, 
the value of E on the contour lines is increasing outwards as desired. 

A bit farther from the curve, i.e. for 1/R < \w\ < r~ a , < a < |, the 
function £ equals asymptotically the one of the corresponding ellipse h^°'(w) — 
rw + ao + r«iw _1 , we denote it by £ (°>. Indeed, remarking that for the area of this 
ellipse we have t =to +0(r 4 ) and that £ (°'(w) — 0{r~ 2a ), we obtain, uniformly 
in w, 



£(w) - £ (0 \w) = -(\h{w)\ 2 - \h (0) (w)\ 2 - \h(l)\ 2 + |h (0) (l)| 
to 



2\ 



9 l' w 

+ -Re/ (hiw-^-h^iw-^h'i^dw 
to J\ 

+ - Re / (h'(w) - ft(°) ' (w))h (a \w- v ) dw 
to J\ 



i(4 0) -*o)^H 



to 
= 0(r l - a ) + 0{r l - 3a ) + 0(r l - 2a ) + 0(r 2 - 2a ) 

-^0 (r -► 0) 

Because £(°\w) > C > for all |w| > i and r > 0, we may choose r so small that 
\£(w) - £(°\w)\ < £(°\w) and so £(w) > for all w E B r - a \ B 1/R . 
It remains the domain, where |ty| > r~ 

h(w) k - (rw) k - 



0(r 2 ). 




And since 



t V{z) = \z\ 2 -t 2 z 2 -t 2 z 2 + o(\z\ 2 ) 



= l -(\z-2t 2 z\ 2 + (l-A\t 2 \ 2 )\z\ 2 )+o(\z\ 2 ) 

for z — > and V > on D \ {0}, we may find a constant C > such that 
toV(z) > C\z\ 2 for all z in the compact domain D. 
Therefore, for r — > 0, 

V(h(w)) = V(rw) + O(l) > —\rw\ 2 + O(l), 

to 

which tends to infinity at least as r~ 2a . 

On the other side, the integral over the logarithm in (15) diverges for r — > only 
as logr. So we have for r small enough E(h(wj) > for all w e h^ 1 (D\D + )\B r - a . 
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This proves now E(z) > for all z E D and therefore, with Corollary 3.5, that 
\x is the equilibrium measure. 

6.4. Shifting the origin. As a little generalization, we consider the case where 
t\ 7^ 0. This corresponds to a shift of the origin. Therefore, we like to define the 
harmonic moments also for a curve which does not encircle the origin. 

Definition. Let 

n 

h(w) = rw + y^ djW^J 

3=0 

parametrize a polynomial curve of degree n. Then the harmonic moments (ife)^^ 
are given by the equation system (9). All other harmonic moments are set to zero. 
Proposition 5.2 tells us that this definition coincides with the previous one if the 
origin is in the interior domain of the curve. 

Corollary 6.7. For any set (ifc)^Li C C with fa\ < \ and tk = for k > n+ 1 and 
any compact domain D C C such that U(z) = \z\ 2 — 2Re^)! =1 t^z k has exactly 
one absolute minimum in the interior of D, there exists a S > so that for all 
< to < S the equilibrium measure jjl for V = j-U on D is given by 

1 

Trio' 

where D + denotes the interior domain of the polynomial curve 7 defined by the 
harmonic moments (ifc)^Lo- 



/»= 3-Xd+A, 



Proof: Let us first shift the origin by ao, such that V(z + ao) has its absolute 
minimum in 0. Thereby, the potential gets the form 



V(z + a ) = U\z\ 2 -2ReJ2t'kA 



V(a ), 



where the t' k are the harmonic moments of the shifted curve 7 — ao. Indeed, the 
coefficients of V and the harmonic moments depend polynomially on the shift ao- 
Because we know them to coincide as long as the origin is inside D + , they do so 
for all ao. 

Applying now Theorem 6.1 for the shifted potential gives the desired result. □ 

7. Discussion 

We have shown that under suitable assumptions on the polynomial p appearing 
in the potential and on the integration range D of the eigenvalues, the asymptotic 
density is uniform with support on a domain uniquely determined by the coefficients 
of the polynomial p. 

It would be interesting to understand what happens at the range of validity 
of our assumptions. If the potential V has more than one minimum in D then 
one should expect for small to to have an equilibrium measure with disconnected 
support, so that a description by a polynomial curve cannot be valid. Also as to 
gets bigger, polynomial curves develop singularities and become non-simple. The 
question is then what happens to the eigenvalues. Finally we note that polynomial 
curves are (real sections of complex) rational curves. Curves of higher genus should 
arise by replacing p by more general holomorphic functions. 
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